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In  Memoriam 
Professor  J.H.  Richmond 
1922  -  1990 


Prof.  J.H.  Richmond  died  suddenly  on  November  30,  1990  of  a  heart 
attack.  This  was  only  two  davs  after  he  completed  the  section  in  ims  report 
on  the  interaction  of  electromagnetic  microwaves  with  nonlinear  devices  and 
nonlinear  media.  Professor  Richmond  has  been  a  member  of  the  staff  at  the 
ElectroScience  Laboratory  since  1952.  He  completed  his  Ph.D.  degree  in 
1955  and  retired  as  a  Professor  from  the  Electrical  Engineering  Department 
of  Ohio  State  University  in  1985.  Prof.  Richmond  was  a  Life  Fellow  in 
the  IEEE  and  was  the  recipient  of  one  of  18  Centennial  awards  from  the 
Antenna  and  Propagation  Society.  He  nevertheless  continued  his  research 
activities  at  the  ElectroScience  Laboratory  and  actually  increased  his  own 
research  activities  during  the  last  few  years  of  his  life. 

Professor  Richmond’s  research  in  the  now  commonly  accepted  moment 
method  represented  a  worldwide  guiding  light  to  users  and  researchers  alike 
in  this  important  technology.  He  represented  an  important  asset  to  the 
ElectroScience  Laboratory  in  many  ways.  As  a  consultant  for  proposed 
theoretical  concepts,  there  was  a  steady  parade  into  his  office  of  senior  staff, 
graduate  students,  and  visitors  to  discuss  their  specific  problems.  The  use 
of  his  computer  codes  was  an  initial  step  for  many  graduate  students.  He 
was  contacted  by  many  external  researchers  who  wanted  to  discuss  some 
phase  of  his  work.  He  will  be  sorely  missed. 

Professor  Richmond’s  published  papers  were  concentrated  in  the  An¬ 
tenna  and  Propagation  Society  of  the  Institute  of  Electrical  and  Electronic 
Engineers  and  Radio  Science.  In  his  later  years,  he  did  not  participate  in 
the  various  electromagnetic  symposia,  so  many  readers  would  not  recognize 
him,  but  all  of  those  involved  in  his  technical  areas  would  be  aware  of  his 
published  research. 
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I.  DIRECTORS  OVERVIEW 


This  report  represents  the  thirteenth  annual  summary  of  The  Ohio  State 
University  Joint  Services  Electronics  Program  (JSEP). 

There  have  been  a  total  27  Ph.D.  and  20  M.Sc  degrees  in  Electrical 
Engineering  obtained  under  partial  JSEP  sponsorship.  There  are  currently 
9  Ph.D.  and  1  M.Sc.  student  being  partially  supported  under  JSEP. 

As  may  be  seen  in  the  Annual  report  Appendix,  14  reprints  have  been 
included  in  the  period  September  1989  to  September  1990.  In  addition,  13 
papers  have  already  been  accepted  for  publication  in  the  coming  year,  an 
additional  12  papers  have  been  submitted,  and  an  additional  12  papers  are 
in  preparation. 

II.  DESCRIPTION  OF  SPECIAL  ACCOM¬ 
PLISHMENTS  AND  TECHNOLOGY  TRAN¬ 
SITION 

i 

The  transfer  of  the  compact  range  and  target  identification  technology  ini¬ 
tiated  under  JSEP  support  for  time  domain  studies  continues  to  make  large 
advances.  Using  other  sources  of  support,  design  for  a  mini  chamber  has 
now  been  generated  and  has  been  constructed.  This  range  is  designed  to 
study  smaller  targets  at  higher  frequencies. 

The  research  has  proven  to  be  of  intense  interest  to  DoD  and  the 
Aerospace  industry  and  our  Compact  Range  Consortium  represents  a  ma¬ 
jor  cross  section  of  the  Aerospace  and  Electronic  Industries,  including  addi¬ 
tional  major  support  from  several  DoD  agencies.  In  fact,  the  total  support 
in  these  experimental  studies  substantially  exceeds  our  JSEP  support.  This 
research  is  truly  guiding  a  major  portion  of  this  technology  in  the  USA  and 
has  been  extremely  important  for  stealth  technology  advances.  However, 
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these  advances  were  only  possible  because  of  the  initial  JSEP  support.  This 
continues  to  be  a  case  where  a  small  investment  of  basic  research  funds  have 
been  leveraged  to  generate  much  larger  support  and  have  achieved  major 
contributions  for  DoD.  This  has  also  lead  to  OSU-ESL  involvement  in  com¬ 
mittees  discussing  Ultra  Wide  Band  radar  systems. 

Our  target  identification  work,  also  partially  funded  at  one  time  under 
JSEP  Time  Domain  Studies,  is  also  being  funded  by  several  other  agencies 
including  ONR  and  continues  to  be  rather  vigorous.  Again,  JSEP  funds 
have  been  leveraged  to  initiate  larger  programs  which  have  been  supported 
continuously  since  JSEP  funding  was  terminated. 

Our  JSEP  research  continues  to  focus  on  Electromagnetic  related  topics. 
There  are  three  major  electromagnetics  areas  that  were  pursued  in  the  past 
year  and  a  closely  related  study  in  Adaptive  Arrays. 

The  goal  of  our  Diffraction  Studies  is  to  not  only  treat  new  diffract¬ 
ing  mechanisms  but  also  to  reduce  their  complexity  so  that  they  can  be 
more  readily  applied  to  DoD  problems.  These  mechanisms  become  exceed¬ 
ingly  important  as  stealth  technology  advances,  i.e.,  as  scattered  fields  are 
reduced  ever  lower.  It  is  the  intent  to  reduce  these  analyses  to  the  Diffrac¬ 
tion  Coefficient  format  so  that  the  solution  of  scattered/radiated  fields  for 
aerospace  vehicles  will  involve  the  use  of  these  coefficients  and  differential 
geometry.  Such  solutions  will  then  be  packaged  in  a.  variety  of  computer 
codes  on  other  projects  as  part  of  a  technology  transfer  mechanism. 

This  activity  has  included  the  generat  ion  and  use  of  simplified  boundary 
conditions  which  are  designated  as  Generalized  Impedance  Boundary  Con¬ 
ditions  (GIBC)  and  Generalized  Resistive  Boundary  Conditions  (GRBC). 
Dr.  Rojas  received  the  best  papej  award  for  ail  IEEE  publications  for  au¬ 
thors  under  35  years  of  age  for  his  effort  in  this  area.  Phis  later  case  is 
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proving  valuable  in  the  design  of  tapered  resistive  cards  for  reducing  scat¬ 
tered  fields  from  edges.  We  are  also  examining  the  effect  of  terminations 
for  open-ended  waveguides  and  this  research  has  resulted  in  a  Ph.D.  disser¬ 
tation.  Our  search  for  an  adequate  corner  diffraction  coefficient  has  been 
rewarding  in  that  an  improved  uniform  corner  diffraction  coefficient  has 
been  generated  and  is  discussed  in  yet  another  dissertation. 

In  the  past,  our  Integral  Equation  Studies  have  focussed  attention  on 
the  analysis  of  penetrable  materials  used  in  conjunction  with  conducting 
surfaces.  These  have  included  various  tactics  to  reduce  the  computation 
time  required,  such  as  making  use  of  special  Green’s  functions  so  that  the 
only  unknowns  are  currents  in  the  materials.  Many  of  the  basic  techniques 
developed  under  .ISEP  support  were  later  incorporated  into  user-oriented 
computer  codes.  These  codes,  such  as  the  “Electromagnetic  Surface  Patch 
Cede  (FSP)”  are  now  widely  used  in  government  and  industry.  Our  current 
research  in  this  area  involves  much  more  general  media  and  we  have  been 
successful  in  treating  chiral  and  non-linear  media  via  the  integral  equations 
approach.  The  research  in  the  chiral  media  scattering  has  resulted  in  several 
papers  and  a  dissertation  is  nearing  completion.  Another  dissertation  by  a 
math  student  is  being  completed  on  the  subject  of  artificial  media. 

The  Hybrid  Approach  represents  novel  analyses  involving  more  than 
one  basic  technique  such  as  was  done  originally  at  the  ESL  by  combining 
diffraction  and  integral  equations  which  was  one  of  the  earlier  such  solu¬ 
tions.  One  of  our  initial  efforts  involved  the  scattering  from  structures 
that  resembled  jet  intakes  and  exhausts.  Several  decades  ago.  these  were 
supposedly  geometries  whose  scattering  properties  would  never  be  treated 
analytically  with  any  degree  of  success.  Our  recent  work  has  been  overcom¬ 
ing  most  of  these  difficulties  as  will  be  seen  in  the  deep  cavities  discussed  in 
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the  appropriate  section.  Both  government  and  industry  are  becoming  the 
primary  supporters  for  this  effort  and  again  JSEP  support  has  been  used  in 
the  initial  stages  of  study  that  have  been  carried  to  t. fie  extent  that,  others 
are  now  providing  the  major  funding.  Research  on  more  shallow  antenna 
cavities  has  continued  under  JSEP  support  and  we  expect  that  this  effort 
will  also  be  of  general  interest  to  many.  We  are  considering  these  same 
techniques  to  evaluate  the  electromagnetics  environment  on  naval  ships. 
This  represents  an  important  topic  for  future  research  and  should  have  a 
substantial  impact  for  the  Navy. 

A  not  her  topic  of  interest  here  is  t  he  elect  romagnetic  properties  of  stripline 
systems.  To  treat  such  devices  rigorously  requires  the  inclusion  of  a  very 
complex  Sommerfeld  integral.  Asymptotic  forms  of  this  integral  have  been 
obtained  that  greatly  simplify  such  analyses.  The  form  of  the  solution  re¬ 
mains  uniformly  valid  even  for  lateral  field  and  source  point  separation  as 
small  as  a  third  of  a  wavelength  and  has  been  extended  to  include  double 
layers.  This  and  future  efforts  should  contribute  substantially  to  the  design 
of  MMIC  st  met  ures. 

t  hese  asymptotic  forms  coupled  to  a  judicious  choice  of  basis  functions 
for  appropriate  choice  of  boundary  conditions  for  the  geometry  (coupler, 
bend,  transformer,  etc.)  not  only  simplify  the  analysis,  but  contribute  sub¬ 
stantially  to  understanding  the  physical  mechanisms  involved.  The  Hybrid 
approach  is  being  used  to  generate  solutions  for  structures  where  neither 
moment  method  or  asymptotics  can  be  expected  to  produce  accurate  an¬ 
swers. 

Our  work  m  array  processing  has  been  concentrated  in  two  areas  this 
year.  First ,  we  have  developed  a  met  hod  for  using  a  multiple-beam  adaptive 
array  (MBAA)  in  both  slotted  and  mishitted  packet  radio  networks,  and 


we  have  evaluated  t  lie  performance  of  single-hop  networks  using  an  MBA  A. 
We  have  shown  that  an  MB  A  A  yields  very  high  throughputs  and  low  delays 
in  such  networks.  Also,  we  have  found  that  the  performance  of  an  unslotted 
network  with  an  MBAA  can  he  almost  as  good  as  that  of  a  slotted  network. 

Second,  we  have  begun  some  new  work  in  angle  estimation.  The  purpose 
of  this  work  is  to  generalize  existing  anele  estimation  techniques  to  make 
it  possible  to  handle  signal  arrival  angles  in  two  spatial  coordinates  and  to 
incorporate  known  signal  modulations  and  unknown  signal  polarizations  in 
the  estimator.  Some  preliminary  results  have  been  obtained. 

Technology  transition  continues  to  take  several  forms  for  our  JSEP  pro¬ 
gram.  First,  of  course,  are  the  students  graduating  in  this  program  who 
carry  the  knowledge  gleaned  in  their  research  programs  to  other  users.  Sec¬ 
ond,  there  are  the  published  papers,  both  oral  and  written,  which  generally 
attract  the  attention  of  other  Dol)  sponsoring  agencies.  Such  agencies  in 
turn  provide  additional  funding  and  in  general  make  use  of  our  JSEP  re¬ 
search  and  extend  it  to  better  their  own  programs. 

Yet  another  method  takes  the  form  of  computer  codes  developed  under 
non-JSEP  sources  that  make  extensive  use  of  JSEP  research.  As  we  have 
noted  previously,  the  results  of  all  of  these  studies  are  of  great  importance 
in  the  analysis  and  control  of  the  radiation  and  scattering  from  complex 
shapes. 

This  continues  to  be  a  major  task  at  The  Ohio  State  University  Elec- 
troScience  Laboratory  (OSU-ESI.)  which  is  funded  bv  a  variety  of  Dol) 
agencies.  A  major  objective  of  the  ESI,  funded  by  other  sources  is  t<>  pro 
vide  a  general  computer  code  (or  codes)  for  the  evaluation  of  the  IK'S  of 
Aerospace  vehicles,  but  a  variety  of  theoretical  analysis  must  be  generated 
before  this  goal  can  become  a  reality.  I  he  OSU-ESI,  continues  to  provide 
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to  DoD  users  a  variety  of  complex  computer  codes  at  the  cost  of  materi¬ 
als  for  radiation  from  antennas  on  aircraft,  reflector  antennas  and  integral 
equation  formulations  based  on  previous  research  activities  to  94  industrial 
organizations  with  DoD  approval  for  use  in  DoD  activities.  In  fact,  last 
year  154  additional  copies  of  these  very  complex  user  friendly  codes  were 
issued.  This  represents  about  a  30%  increase  over  previous  years  in  spite 
of  cutbacks  of  DoD  funding. 
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III.  DIFFRACTION  STUDIES 

Researchers: 

P.H.  Pathak,  Associate  Professor 
R.C.  Chou,  Assistant  Professor 
R.  Rojas,  Senior  Research  Assoc. 

R.J.  Burkholder,  Post-Doctoral  Researcher 

G.  Zogbi,  Graduate  Research  Assoc. 

H. C.  Ly,  Graduate  Research  Assoc. 

L.M.  Chou,  Graduate  Research  Assoc. 

1.  Introduction 

The  topics  investigated  during  the  past  year  consisted  of  the  electromag¬ 
netic  (EM)  diffraction  by  a  corner  in  a  plane  metallic  angular  sector,  the 
EM  diffraction  by  discontinuities  in  planar,  dielectric/magnetic  material 
boundaries,  the  high  frequency  EM  field  representations  in  terms  of  beams, 
or  rays,  for  electrically  large,  closed  guiding  structures,  and  the  represen¬ 
tation  of  EM  caustic  fields  in  a  form  suitable  for  analyzing  the  reflection 
and  diffraction  of  fields  generated  by  smooth  caustics  near  the  points  of 
reflection  and  diffraction.  A  substantial  progress  has  been  made  on  these 
topics  as  described  below. 

2.  EM  Diffraction  by  a  Corner  in  a  Plane  Angular 
Sector 

A  uniform  geometrical  theory  of  diffraction  (UTD)  solution  has  been  ob¬ 
tained  for  describing  the  diffraction  by  a  corner  in  a  planar  perfectly- 
conducting  surface.  This  UTD  corner  diffraction  solution  constitutes  an 
important  contribution  to  the  UTD  area  in  that  a  UTD  solution  for  the 
diffraction  by  such  a  fundamental  geometry  has  not  been  developed  previ- 
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ously,  except  for  one  which  was  conjectured  heuristically  [1]  and  which  is  not 
as  versatile  and  accurate  as  the  present  solution.  This  solution  is  uniform 
in  the  sense  that  it  keeps  the  total  high  frequency  field  continuous  across 
the  corner  induced  edge  diffracted  ray  shadow  boundary  cones  where  the 
edge  diffracted  fields  become  discontinuous.  Continuity  of  the  fields  is  es¬ 
tablished  here  in  the  same  manner  as  that  done  by  the  UTD  edge  diffracted 
field  which  makes  the  total  wedge  diffracted  field  continuous  across  the  edge 
induced  geometrical  optics  incident  and  reflected  ray  shadow  boundaries. 
The  discontinuities  in  the  geometrical  optics  incident  and  reflected  fields  at 
the  incident  and  reflected  ray  shadow  boundaries,  that  are  formed  by  the 
presence  of  edges  of  the  plane  angular  sector,  are  compensated  by  the  fields 
diffracted  by  these  edges;  whereas,  the  discontinuities  in  the  edge  diffracted 
fields  past  the  corner  are  compensated  by  the  corner  diffracted  field.  The 
special  function,  or  the  transition  function,  which  allows  the  edge  diffracted 
fields  to  compensate  for  the  discontinuities  in  the  geometrical  optics  field 
involves  the  Fresnel  integral  as  is  well  known  in  the  UTD  for  edges.  A 
different  transition  function  is  necessary  which  allows  a  proper  compensa¬ 
tion  of  the  discontinuities  of  the  edge  diffracted  fields  past  the  corner;  this 
transition  function  has  been  found  in  the  present  work.  Furthermore,  the 
properties  of  this  uniform  corner  diffraction  transition  function  have  been 
studied  carefully,  both  for  small  arguments  (i.e.  near  the  corner  induced 
edge  diffraction  shadow  boundaries)  and  large  arguments  (i.e.  far  from  the 
corner  induced  edge  diffraction  shadow  boundaries).  Most  importantly,  a 
very  efficient  method  and  an  algorithm  based  on  if  have  been  developed 
to  evaluate  the  UTD  corner  transition  function;  Ibis  slep  is  crucial  in  f  he 
use  of  this  UTD  corner  diffraction  solution  to  solve  practical  diffraction 
problems.  While  this  UTD  corner  solution  is  based  on  the  best  available 
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estimate  for  the  asymptotic  approximation  for  the  current  induced  on  the 
plane  angular  sector,  it  provides  results  that  are  in  very  good  agreement 
with  other  rigorous  numerical  approaches  for  treating  the  EM  scattering 
by  finite  plates  with  corners.  The  present  solution  is  valid  for  aspects  not 
close  to  grazing  where  the  launching  and  diffraction  of  edge  waves  produced 
by  the  corner  become  important;  this  is  due  to  the  fact  that  the  present 
asymptotic  approximation  for  the  current  on  the  plane  angular  sector  is 
unable  to  model  the  effects  of  edge  waves  resulting  from  the  presence  of  the 
corner.  The  inclusion  of  the  latter  effect  is  a  topic  of  future  research;  some 
initial  attempts  at  predicting  these  effects  have  also  been  made  recently  at 
our  laboratory. 

3.  Diffraction  by  Discontinuities  in  Planar  Material 
Boundaries 

As  reported  in  the  1989  annual  report,  the  study  of  high  frequency  scatter¬ 
ing  by  non-conducting  and  penetrable  surfaces  can  be  somewhat  simplified 
by  developing  equivalent  boundary  conditions  for  characterizing  such  sur¬ 
faces.  These  boundary  conditions  are  referred  to  as  Generalized  Impedance 
(for  opaque  bodies)  and  Resistive  (for  transparent)  Boundary  Conditions 
(GIBC  and  GRBC).  Once  these  boundary  conditions  are  developed,  it  is 
no  longer  necessary  to  calculate  the  fields  inside  the  material  and  they  are 
also  very  suitable  in  the  analysis  of  scattering  problems  by  means  of  various 
functional  analytical  techniques.  One  disadvantage  is  the  fact  that  the  so¬ 
lutions  developed  with  these  boundary  conditions  are  in  general  not  unique. 
A  lot  of  effort  has  been  directed  to  develop  additional  conditions  such  as 
edge  and  junction  conditions  which  are  more  general  than  the  usual  edge 
condition.  In  particular,  a  junction  condition  has  been  developed  at  the 
planar  junction  of  two  thin  dielectric  half-planes.  Note  that  the  application 
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of  this  junction  condition  yields  a  unique  solution. 

A  more  general  type  of  coating  can  be  built  with  chiral  material.  This 
type  of  material  has  been  known  to  exist  at  optical  frequencies,  however, 
it  appears  that  it  can  be  built  at  microwave  frequencies.  A  chiral  material 
can  be  described  bj  three  parameters,  namely,  its  permittivity,  permeabil¬ 
ity  and  an  additional  parameter,  referred  to  as  the  chiral  admittance.  Thus, 
there  is  one  more  degree  of  freedom  which  can  be  used  to  control  the  scat¬ 
tering  properties  of  a  chiral  body.  To  study  the  EM  diffraction  by  bodies 

* 

constructed  with  chiral  material,  GIBC/GRBC’s  have  been  developed  for 
planar  chiral  slabs,  which  can  be  grounded  by  a  perfect  electric  conductor. 
As  in  the  case  of  dielectric  materials,  these  equivalent  boundary  conditions 
will  facilitate  the  analysis  of  diffraction  problems. 

It  is  customary  to  compare  newly  developed  diffraction  solutions  with 
other  independent  solutions  such  as  those  based  on  the  moment  method 
technique.  For  these  purposes,  integral  equations  have  been  developed  for 
two-  and  three-dimensional  chiral  bodies.  These  integral  equations  are  more 
efficient  than  the  ones  developed  by  a  straight  forward  application  of  the 
volume  equivalent  current  method. 

4.  Ray  and  Beam  Representations  for  Fields  Coupled 
into  Large  Open  Cavities 

A  Gaussian  beam  shooting  approach  was  developed  for  analyzing  the  EM 
fields  coupled  into  an  electrically  large  open  cavity  of  relatively  arbitrary 
shape,  when  the  open  end  is  directly  illuminated  bv  an  external  source. 
The  fields  coupled  into  the  cavity  via  the  open  end  are  described  by  beams 
which  are  launched  from  an  array  of  points  in  the  aperture  at  the  open 
end.  This  array  of  launch  points  corresponds  to  the  phase  centers  of  the 
subapertures  into  which  the  original  aperture  is  subdivided.  A  sufficiently 
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dense  set  of  identical  beams  is  launched  radially  from  each  phase,  center 
and  with  constant  interbeam  angular  spacing.  The  initial  beam  launching 
amplitudes  are  determined  by  matching  the  beams  from  each  subaperture 
to  the  far  field  pattern  of  that  subaperture  with  the  cavity  walls  removed. 
Furthermore,  once  launched,  the  beams  are  tracked  through  the  interior 
cavity  region  using  an  axial  tracking  approximation.  In  this  approximation, 
the  beams  are  tracked  like  real  rays  along  their  axes  using  the  rules  of 
beam  optics.  To  use  this  approximation,  it  is  necessary  to  choose  highly 
focussed  (or  wide  waisted)  beams.  However,  such  beams  undergo  distortion 
upon  each  reflection  at  the  interior  walls.  Thus,  this  real  ray  axial  beam 
tracking  approximation  has  been  found  to  work  well  provided  the  length 
of  the  waveguide  cavity  is  not  much  larger  than  four  or  five  times  the 
maximum  aperture  dimension  at  the  open  front  end.  One  could  overcome 
this  limitation  by  tracking  beams  via  complex  rays;  however,  complex  ray 
tracking  appears  at  this  time  to  be  a  very  cumbersome  and  time  consuming 
procedure.  A  significant  advantage  of  the  beam  approach  for  representing 
the  interior  cavity  field  is  that  the  beams  need  to  be  tracked  only  once  into 
the  interior  cavity  from  the  array  of  points  in  the  aperture  at  the  open  end 
which  is  illuminated.  The  beam  paths  then  remain  unchanged  with  change 
of  illumination  or  excitation;  only  the  initial  beam  amplitudes  change  with 
changes  in  illumination  of  the  aperture  at  the  open  end  which  excites  the 
cavity.  Furthermore,  the  beams  launched  into  the  cavity  automatically  and 
implicitly  include  the  effects  of  diffraction  at  the  boundary  (edge)  of  the 
aperture. 

An  alternative  to  tracking  beams  via  complex  rays  to  overcome  the 
effects  of  beam  distortion  indicated  earlier  is  to  launch  a  dense  grid  of  ordi¬ 
nary  ray  tubes  radially  from  the  subapertures  instead  of  launching  beams. 
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Once  launched,  the  far  thinner  ray  tubes  are  then  tracked  via  reflections 
through  the  cavity  via  the  rules  of  ray  optics.  This  procedure  is  referred 
to  as  the  generalized  ray  expansion  (GRE)  which  has  been  developed  and 
shown  to  work  quite  well  for  several  representative  cases  that  have  been 
studied  recently.  This  recently  developed  GRE  method  maintains  the  use¬ 
ful  features  of  the  beam  expansion  in  that  the  ray  paths  in  the  GRE  need 
to  be  tracked  only  once  because  they  do  not  change  with  the  illumination, 
and  the  GRE  also  implicitly  contains  effects  of  the  fields  which  couple  into 
the  cavity  via  diffraction  by  the  edges  of  the  aperture  at  the  open  end.  It  is 
also  noted,  however,  that  a  lot  more  rays  need  to  be  tracked  with  GRE  than 
one  needs  to  track  beams  in  the  beam  shooting  approach.  It  is  noted  that 
the  conventional  Geometrical  Optics  ray  shooting  schemes  require  a  new 
set  of  incident  rays  entering  the  cavity  to  be  tracked  each  time  the  incident 
angle  changes,  and  they  must  be  modified  to  explicitly  include  the  effects 
of  fields  coupled  into  the  cavity  via  diffraction  at  the  open  end  by  tracking 
within  the  cavity  a  whole  additional  set  of  rays  diffracted  by  the  aperture 
edges,  the  GRE  on  the  other  hand  does  not  have  these  deficiencies. 

As  a  part  of  a  future  effort  on  this  topic,  some  other  ray  and  beam 
launching  schemes  will  be  studied  which  would  be  more  easily  and  conve¬ 
niently  applicable  to  an  aperture  of  relatively  arbitrary  shape  at  the  open 
end.  Additional  refinements  to  the  GRE  will  also  be  investigated  to  improve 
the  efficiency  of  the  GRE  procedure. 

5.  Caustic  Field  Analysis 

The  object  here  is  to  find  the  field  scattered  from  an  object  when  the  il¬ 
luminating  field  is  characterized  by  a  smooth  caustic  of  rays  in  which  the 
caustic  is  very  close  to  the  scattering  surface.  If  the  caustic  is  far  from 
the  scatterer,  then  the  field  away  from  the  caustic  transition  region  may  be 
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represented  as  a  ray  optical  field  whose  scattering  can  thus  be  calculated 
using  ordinary  ray  methods.  Previous  work  on  our  JSEP  yielded  an  asymp¬ 
totic  high  frequency  solution  for  the  field  near  a  smooth  caustic  which  is 
uniformly  valid  across  the  caustic  and  which  reduces  to  the  correct  real  and 
complex  ray  optical  solution  in  the  lit  and  dark  regions,  respectively,  out¬ 
side  the  caustic  transition  region.  However,  it  is  not  possible  to  directly  use 
that  uniform  expresssion  to  perform  an  asymptotic  high  frequency  calcula¬ 
tion  of  the  subsequent  reflection  or  diffraction  of  this  caustic  field,  because 
it  is  not  in  the  form  of  some  type  of  superposition  of  simpler  ray  optical 
fields.  What  is  therefore  needed  is  to  express  the  field  near  a  caustic  as  a 
superposition  of  simple  ray  optical  constituents  which  can  then  be  tracked 
upon  reflection  and  diffraction  via  ray  optics.  With  the  above  objective 
in  view,  a  plane  wave  (spectral)  representation  has  been  obtained  for  the 
caustic  field  generated  by  a  concave  reflector.  It  is  next  desired  to  reflect 
or  diffract  this  field  by  a  second  scattering  object  to  ascertain  the  field 
reflected  or  diffracted  by  the  latter  when  illuminated  by  this  caustic  field. 
This  part  of  the  study  is  now  in  progress.  It  is  noted  that  a  complication 
can  arise  because  the  fields  reflected  or  diffracted  from  the  second  scatterer 
could  also  exhibit  a  smooth  caustic;  this  situation  will  also  be  incorporated 
in  the  study. 
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1.  Introduction 

This  section  will  summarize  our  work  in  integral  equation  studies  from  July 
1,  1989  to  the  present.  In  overview,  our  recent  research  has  centered  on 
integral  equation  and  method  of  moments  (MM)  solutions  for  unusual  or 
novel  media.  In  particular,  we  are  developing  MM  solutions  for  scattering 
by  chiral  bodies,  an  MM  solution  for  the  chiral  microstrip  transmission  line, 
and  also  MM  solutions  for  the  analysis  of  artificial  media  and  nonlinear 
media. 

a.  Chiral  Media 

A  chiral  media  has  constitutive  relationships  of  the  form 

D  tE  ji  B 
H  -  — B  ji E 

where  £  is  the  chiral  admittance.  In  a  regular  or  achiral  media,  £  —  0.  The 
non-zero  (  in  a  chiral  media  results  in  an  additional  coupling  between  I  lie 
electric  and  magnetic  fields,  and  causes  the  polarization  of  a  plane  wave  to 
rotate  as  it  propagates  through  the  chiral  media. 

One  component  of  our  research  in  chiral  media  is  directed  toward  the 
development  of  an  integral  equation  and  MM  solution  for  the  scattering 
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from  a  chiral  cylinder  of  arbitrary  cross  section  shape.  In  order  to  treat 
inhomogeneous  chiral  cylinders  we  employ  the  volume  current  formulation 
in  which  the  chiral  cylinder  is  represented  by  equivalent  electric  and  mag¬ 
netic  volume  polarization  currents  (J,M).  To  this  end,  we  have  developed 
a  volume  equivalence  theorem  for  chiral  media,  and  have  used  this  theorem 
to  formulate  a  pair  of  coupled  vector  integral  equations  for  the  equivalent 
electric  and  magnetic  currents  [1,2].  These  coupled  vector  equations,  which 
are  equivalent  to  six  coupled  scalar  equations,  are  solved  via  a  pulse  basis 
and  point  matching  MM  solution. 

In  order  to  verify  the  accuracy  of  the  MM  solution,  we  developed  an 
(new)  exact  eigenfunction  solution  for  scattering  by  a  multilayered  circular 
chiral  cylinder  [3] .  For  example,  Figure  1  shows  a  comparison  of  the  MM 
and  eigenfunction  solution  for  TE  scattering  by  a  circular  chiral  cylinder. 
Numerical  results  generated  by  the  MM  solution  are  in  excellent  agreement 
with  the  exact  eigenfunction  solution.  Note  also  that  the  chirality  of  the 
media  can  significantly  effect  the  scattering,  with  a  major  efTect  being  a  ro¬ 
tation  of  the  polarization  of  the  scattered  wave.  Recently  we  have  modified 
the  above  MM  solution  to  include  the  presence  of  a  perfectly  conducting 
half  plane  [1]. 

Image  theory  is  a  simple  technique  for  analyzing  bodies  over  a  ground 
plane.  We  have  extended  conventional  image  theory  to  include  chiral  bodies 
[5],  As  illustrated  in  Figure  2,  the  image  of  a  chiral  body  obeys  conventional 
image  theory  except  for  a  minus  sign  in  the  chiral  admittance. 

As  illustrated  in  Figure  3,  our  most  recent  work  in  chiral  media  involves 
a  microstrip  transmission  line  on  a  chiral  substrate.  I  bis  work  required 
the  development  of  tfie  spectral  domain  Green's  function  for  a  grounded 
chiral  slab.  This  chiral  slab  Green's  function  is  then  incorporated  into  a 
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Figure  1:  A  comparison  of  the  MM  and  eigenfunction  solutions  for  TE 
scattering  by  a  chiral  cylinder. 
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CHIRAL  IMAGE  THEORY 


Figure  2:  Image  theory  for  a  chiral  body. 
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Figure  3:  A  microstrip  transmission  line  on  a  chiral  substrate. 
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ARTIFICIAL  DIELECTRIC/FERRITE  MEDIA 


Figure  4:  An  artificial  medium  is  modeled  by  a  periodic  array  of  small 
scatterers. 

MM  solution  for  the  current  distribution  J(x)  and  propagation  constant  f3 
of  a  microstrip  transmission  line  on  the  chiral  substrate. 

b.  Artificial  Media 

As  illustrated  in  Figure  4,  an  artificial  medium  is  created  by  suspending  a 
large  number  of  small  scatterers,  such  as  spheres,  discs,  or  dipoles,  in  some 
host  or  background  medium.  A  field  propagating  in  1  his  artificial  medium 
will  induce  electric  and/or  magnetic  dipole  moments  in  the  scatterers.  The 
result  is  that  the  artificial  medium  appears  to  be  a  dielectric  and/or  a 
ferrite  medium.  An  artificial  chiral  media  can  be  created  from  a  3D  infinite 
periodic  array  of  helices. 
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Our  research  is  directed  toward  the  use  of  integral  equation  and  MM 
techniques  for  the  analysis  of  artificial  media.  Essentially,  one  uses  the 
periodic  MM  to  obtain  the  impedance  matrix  for  an  infinite  periodic  array 
of  scatterers  filling  all  space.  By  setting  the  determinant  of  this  matrix 
to  zero  one  can  determine  the  normal  modes  of  the  system  and  in  turn 
deduce  the  equivalent  permittivity  and  permeability  of  the  artificial  media. 
At  present  we  have  applied  this  method  to  the  2D  problem  of  an  infinite 
periodic  array  of  thin  dielectric  cylinders,  and  the  3D  problem  of  an  infinite 
periodic  array  of  short  dipoles.  The  results  of  this  research  will  be  described 
in  a  PhD  dissertation  now  in  preparation  [7]. 

c.  Nonlinear  Media 

We  are  interested  in  the  interaction  of  electromagnetic  microwaves  with 
nonlinear  devices  and  nonlinear  media.  Although  such  interactions  have 
potentially  useful  applications,  very  little  effort  is  currently  devoted  to  this 
area.  In  the  past  37  years,  the  nonlinear  area  has  been  represented  in  the 
IEEE  Transactions  on  Antennas  and  Propagation  only  by  a  few  papers  on 
nonlinear  plasmas,  EMP  transmission  through  a  ferromagnetic  shield,  and 
scattering  by  a  single  dipole  with  nonlinear  load.  On  the  other  hand,  the 
areas  of  nonlinear  optics  and  nonlinear  circuits  show  considerable  activity. 

A  section  on  nonlinear  media  was  included  in  The  1989  URSI  Interna¬ 
tional  Symposium  on  Electromagnetic  Theory.  This  section  consists  of  four 
papers  from  Russia,  one  paper  each  from  Greece,  Italy,  Poland,  Scotland 
and  Sweden,  but  none  from  the  U.S.A.  Out  of  these  nine  papers,  three  are 
concerned  specifically  with  nonlinear  optics. 

We  have  investigated  plane- wave  scattering  by  a  periodic  planar  array 
of  parallel  wire  dipoles.  As  illustrated  in  Figure  5  the  single  row  of  dipoles 
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Figure  5:  One-row  periodic  planar  array  of  dipoles  with  nonlinear  loads. 


forms  a  broadside  array,  and  each  dipole  has  a  nonlinear  load  at  the  center. 
The  load  may  be  series-opposing  diodes  to  represent  a  nonlinear  capacitor, 
or  parallel-opposing  diodes  to  represent  a  nonlinear  resistor.  Several  time- 
harmonic  plane  waves  may  simultaneously  be  incident  on  the  array,  but 
each  wave  has  the  same  direction  of  propagation  and  the  frequencies  are 
harmonically  related.  In  the  steady-state  condition  the  total  field  is  thus 
periodic  in  space  and  time,  and  Floquet’s  theorem  is  applicable. 

Even  though  the  incident  wave  contains  only  one  or  several  harmonic 
frequencies,  the  scattered  field  will  contain  an  infinite  series  of  harmonics. 
This  series  may  be  truncated,  however,  and  we  wish  to  solve  numerically 
for  the  complex  amplitudes  of  the  various  harmonic  plane  waves  scattered 
(or  reflected)  in  the  specular  direction.  One  may  also  be  interested  in  the 
amplitudes  of  some  of  the  grating  lobes.  It  has  been  observed  that  new- 
grating  lobes  may  appear  when  linear  loads  are  replaced  with  nonlinear 
loads. 

The  solution  of  this  problem  is  simplified  by  splitting  it  into  two  parts. 


In  Tart  I  the  nonlinear  loads  are  removed,  and  we  apply  the  moment  method 
to  solve  the  scattering  problem  where  a  time-harmonic  plane  wave  has 
oblique  incidence  on  a  periodic  planar  array  of  dipoles  with  short-circuited 
terminals.  This  procedure  is  carried  out  for  each  harmonic  wave  incident 
on  the  array.  For  the  nth  harmonic,  one  calculates  and  stores  the  short- 
circuit  terminal  current  or  the  open-circuit  received  voltage  Tnr  induced  in 
the  dipoles. 

Part  I  also  includes  the  solution  of  a  phased-array  transmitting  problem. 
Again  the  nonlinear  loads  are  removed,  and  a  one- volt  generator  is  placed  at 
the  terminals  of  each  dipole.  The  moment  method  is  then  applied  to  solve 
this  time-harmonic  transmitting  problem.  This  procedure  is  carried  out  for 
each  harmonic  frequency  deemed  to  be  significant  in  the  representation  of 
the  total  scattered  field.  For  the  nth  harmonic,  one  calculates  and  stores 
the  active  antenna  impedance  Zn. 

Figure  6  illustrates  the  Thevenin  equivalent  circuit  for  a  dipole  array 
with  nonlinear  loads.  The  open-circuit  received  voltage  vr(t)  and  the  load 
voltage  vi(t)  are  expanded  in  Fourier  series  as  follows: 

vr(t)  =  Real  ^2  exp(jnu;0<)  (1) 

n 

v,(t)  -  Real^TV^expljnvJ)  (2) 

n 

Of  course  the  active  antenna  impedance  Za( u;)  is  frequency  dependent,  and 
Zn  denotes  its  value  at  the  nih  harmonic  frequency.  From  the  basic  circuit 
laws  as  applied  to  Figure  6,  the  load  current  is  given  by 

*(0  =  RealJ^  “^4,— exp(jnw„<)  ~  ?,(«)  (3) 

n 

Now  the  load  current  i(t)  and  the  load  voltage  \){t)  must  be  related 
by  the  characteristics  of  the  nonlinear  load.  If  each  dipole  is  loaded  with 
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Figure  6:  Thevenin  equivalent  circuit  for  dipole  array  with  nonlinear  loads. 

matched  parallel-opposing  diodes,  the  current-versus- volt  age  law  may  be 
represented  as  follows: 

i(t)  =  £  cnlvt(t)lnsgn(t)  =  i2(t)  (4) 

n=l 

where  sgn(t)  =  1  when  t  is  positive,  and  -1  when  t  is  negative.  The  analysis 
may  be  extended  to  take  account  of  the  load  inductance  and  capacitance, 
but  this  detail  will  not  be  covered  here. 

At  the  outset,  one  is  given  the  angle  of  incidence,  the  frequency,  and 
the  complex  amplitude  of  each  incident  plane  wave.  The  impedances  Zn 
and  the  received  voltages  V*  are  then  calculated  as  Part  1  of  the  analysis. 
The  object  of  Part  II  is  to  determine  numerically  the  load  voltage  V„.  This 
is  accomplished  as  follows. 

From  Equations  (3)  and  (4),  the  currents  i\{t)  and  ?’ 2 ( < )  must  be  equal  to 
each  other  and  to  i(t).  In  a  least-square  solution,  one  adjusts  the  voltages 
to  minimize  the  quantity  F  which  is  defined  as  follows: 

F  —  f  MO  -  i2(0)2  dt  (5) 

Jo 


26 


where  T  =  27t/u i0  denotes  the  basic  period  of  the  waveform.  We  have  found 
the  gradient  search  technique  to  be  accurate  and  efficient  for  minimizing 
F.  After  the  load  voltages  and  load  currents  In  have  been  determined  in 
this  manner,  it  is  a  simple  matter  to  calculate  the  amplitude  of  the  reflected 
plane  wave  and  the  grating  lobes  scattered  by  the  periodic  structure. 

Figure  7  illustrates  the  backscatter  echo  width  for  a  periodic  array  of 
short-circuited  dipoles  with  broadside  incidence,  as  calculated  with  our  new 
computer  program  which  was  developed  as  part  of  this  study.  For  compar¬ 
ison  Figure  7  also  shows  the  exact  backscatter  echo  width  for  a  perfectly 
conducting  strip  with  infinite  length,  and  the  close  agreement  between  the 
two  results  tends  to  verify  the  accuracy  of  the  calculated  data.  In  Figure  7 
the  wire  dipoles  have  perfect  conductivity,  a  radius  of  0.005  wavelength 
and  a  spacing  of  0.05  wavelength,  and  the  width  of  the  strip  is  equal  to 
the  dipole  length.  Of  course,  when  the  spacing  is  increased,  the  backscat¬ 
ter  characteristics  of  the  dipole  array  will  no  longer  resemble  those  of  the 
conducting  strip. 

For  the  next  example,  we  selected  a  nonlinear  load  consisting  of  a  pair  of 
Schotlky  diodes  (type  HSCH-3486)  connected  in  parallel-opposing.  Thus, 
one  diode  conducts  when  the  load  voltage  is  positive,  and  the  other  conducts 
when  the  load  voltage  is  negative.  Figure  8  illustrates  the  typical  current- 
versus-voltage  characteristics  for  the  HSCH-3486  diode,  together  with  a 
simple  three-term  approximation. 

Figure  9  illustrates  the  load  current  waveforms  i\(t)  and  ?2(f)  defined 
in  Equations  (3)  and  (4).  The  close  agreement  between  the  two  calculated 
waveforms  tends  to  verify  the  accuracy  of  the  solution.  In  this  example  two 
harmonic  plane  waves  have  broadside  incidence  simultaneously  on  the  one- 
r°w  periodic  array  of  dipoles.  One  incident  wave  has  0.3  volts  per  meter  and 
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Figure  7:  Broadside  backscatter  versus  length  for  dipole  array  and  infinite 
strip. 
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Figure  8:  Current  versus  voltage  for  Schottkv  diode  HSCH-3486. 


29 


a  frequency  of  600  MHz,  while  the  other  has  0.15  v/m  at  900  MHz.  At  the 
fundamental  frequency  (300  MHz),  the  wire  radius  is  0.0002A,  the  dipole 
length  is  0.25A,  and  the  dipole  spacing  is  0.04A.  The  free-space  wavelength 
is  denoted  by  A,  and  the  dipoles  are  constructed  of  brass  with  a  conductivity 
of  10.4  MS/m.  Ten  harmonic  frequencies  (300  MHz  to  3  GHz)  were  retained 
in  calculating  the  load  current  and  the  load  voltage.  The  scattered  field 
includes  a  significant  300  MHz  component,  which  may  be  considered  a 
subharmonic  since  this  frequency  is  not  present  in  the  incident  wave. 

To  obtain  accurate  results,  the  sinusoidal-Galerkin  moment  method  was 
employed  together  with  the  magnetic- frill  model  of  the  terminal  region  of 
each  loaded  dipole.  At  the  highest  frequency,  the  length  of  each  dipole  is 
2.5A  and  the  moment  method  involves  23  equations  with  23  unknowns. 

In  this  example  the  gradient-search  techniques  optimized  twenty  param¬ 
eters,  namely  the  real  and  imaginary  parts  of  ten  harmonic  load  voltages. 
In  this  way  it  minimized  the  error  between  the  current  waveforms  i\(t)  and 
*j(0- 

In  the  next  period,  data  will  be  compiled  to  illustrate  some  of  the  in¬ 
teresting  results  obtainable  with  this  configuration,  and  a  paper  will  be 
submitted  for  publication.  In  addition,  an  attempt  will  be  made  to  inves¬ 
tigate  the  solutions  which  may  propagate  on  the  one-row  periodic  array  of 
dipoles. 

2.  Other  Research 

The  results  of  some  of  our  past  work  were  published  during  (his  reporting 
period.  This  work  includes: 

1.  an  eigenfunction  solution  for  scattering  by  a  material  coated  parabolic 
cylinder  [6], 
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2.  an  invited  review  article  on  MM  modeling  of  complex  bodies  [8], 

3.  an  invited  tutorial  review  article  on  the  MM  has  been  prepared  [9]. 
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1.  Introduction 

During  the  past  year,  the  topics  investigated  were  in  the  areas  of  passive 
MMIC  elements  and  antennas,  in  the  hybrid  formulation  of  EM  scatter¬ 
ing  from  and  penetration  into  electrically  large  open  cavities  containing  a 
complex  interior  obstacle,  in  the  hybrid  formulation  of  EM  scattering  by 
shallow  antenna  cavities,  and  in  the  hybrid  combination  of  high  and  low 
frequency  methods  for  predicting  the  EM  scattering  by  complex  objects. 
Considerable  progress  on  these  topics  has  been  made  as  described  below. 

a.  Hybrid  Analysis  of  Microstrip  and  Related  Configurations 

In  order  to  develop  an  efficient  hybrid  analysis  of  passive  MMIC  elements 
and  antennas,  it  is  important  to  study  the  guiding/radiation  properties  of 
various  MMIC  transmission  line  configurations.  A  hybrid  analysis  leading 
to  a  characterization  of  the  MMIC  components  and  antennas  can  then  be 
formulated  in  terms  of  an  integral  equation  which  can  be  solved  via  the 
moment  method  (MM),  such  an  MM  employs  a  hybrid  set  of  basis  func¬ 
tions  for  describing  the  unknown  currents  on  the  MMIC  structures  which 
typically  reside  on  a  grounded  substrate.  In  particular,  a  hybrid  basis  set 
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includes  the  conventional  subsectional  type  in  the  discontinuity  region  of 
the  MMIC  elements  and  the  traveling  wave  guiding/radiating  modes  on  the 
rest  of  the  appropriate  regions  (typically  corresponding  to  the  transmission 
line  regions).  Furthermore,  in  open  region  problems,  the  conventional  form 
of  the  Green’s  function  in  the  kernel  of  the  integral  equation  becomes  highly 
inefficient  for  numerical  computations;  therefore,  an  asymptotic  form  of  the 
planar  Green’s  function  with  the  asymptotic  large  parameter  being  the  lat¬ 
eral  separation  of  the  source  and  field  points  (along  the  substrate)  has  been 
developed,  such  that  it  remains  valid  uniformly  even  for  lateral  field  and 
source  from  separations  as  small  as  a  third  of  the  free  space  wavelength. 
Such  an  asymptotic  Green’s  function  was  previously  developed  on  our  JSEP 
for  the  single  layer  substrate.  More  recently,  it  has  been  extended  to  a  dou¬ 
ble  layer  (superstrate-substrate  case)  and  is  currently  being  extended  to  a 
triple  layer.  The  merits  of  extending  the  asymptotic  analysis  to  arrive  at  a 
Green’s  function  which  is  valid  for  N  number  of  layers  with  N  >  3  are  not 
clear  at  this  time;  this  aspect  is  under  study.  In  addition  to  open  region 
structures,  closed  region  structures,  are  also  currently  under  investigation, 
and  the  utility  as  well  as  the  limitations  of  employing  closed  region  configu¬ 
rations  to  provide  information  on  open  region  ones  will  also  be  subsequently 
investigated;  this  connection  between  the  open  and  closed  regions  has  not 
been  adequately  studied  in  the  past  except  for  the  fundamental  trapped 
mode  on  some  open  guiding  structures.  For  modes  that  leak  energy  from 
the  transmission  line  structure  (e.g.  microstrip  line,  co-planar  line,  slot  line, 
etc.)  the  ability  to  use  a  closed  region  structure  does  not  become  immedi¬ 
ately  apparent.  All  of  the  above  studies  serve  as  building  blocks  wherefrom 
efficient  analysis  and  design  of  a  variety  of  MMIC  structures  and  antennas 
can  be  performed. 
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b.  EM  Scattering  from  and  Penetration  into  Electrically  Large 
Open  Cavities  containing  a  Complex  Interior  Obstacle 

The  EM  scattering  from  large,  open  cavity  structures  containing  a  complex 
obstacle  is  of  great  interest  in  EM  scattered  field  as  well  as  EMC  applica¬ 
tions.  The  cavity  shape  can  be  relatively  arbitrary.  This  coupled  with  the 
fact  that  the  interior  obstacle  can  be  complex  makes  the  study  of  this  scat¬ 
tering  configuration  a  rather  difficult  task.  However,  it  is  found  that  one 
can  break  up  the  analysis  of  such  a  general  cavity-obstacle  configuration 
into  two  distinct  physical  parts.  This  is  done  by  conveniently  choosing  a 
plane  St  near  the  obstacle.  This  plane  divides  the  region  inside  the  cavity 
into  one  from  the  open  front  end  of  the  cavity  to  the  plane  St,  and  into 
another  from  St  to  the  back  end  of  the  cavity  which  contains  the  complex 
obstacle.  The  two  parts  are  analyzed  separately;  i.e.,  as  two  physically 
different  situations  distinct  from  each  other.  In  the  first  part,  the  fields 
are  tracked  within  the  first  section  of  the  cavity  from  the  open  front  end, 
which  is  directly  illuminated  by  an  external  source,  to  the  plane  St',  these 
fields  are  then  allowed  to  radiate  out  of  St  using  a  Kirchhoff  approxima¬ 
tion  that  is  accurate  for  electrically  large  cavities.  In  the  second  part,  a 
unit  amplitude  plane  wave  EM  scattering  characterization  of  the  obsta¬ 
cle  is  obtained  for  the  second  section  from. Sr  to  the  back  of  the  cavity; 
this  section  contains  the  obstacle.  Such  a  plane  wave  scattering  character¬ 
ization,  where  the  response  to  plane  waves  coming  from  a  set  of  different 
directions,  can  be  developed  either  by  analytical,  or  by  numerical  methods 
(e.g.  moment  method  solution  of  the  integral  equation  for  scattering  by 
the  second  section,  or  by  some  other  numerical  procedure  to  achieve  the 
same)  if  the  complex  obstacle  can  be  modeled  adequately  in  this  numerical 
approach,  or,  in  principle,  by  measurements  if  the  obstacle  is  too  complex 
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to  be  treated  using  numerical  approaches.  The  information  based  on  the 
plane  wave  scattering  characterization  of  the  obstacle  part  of  the  cavity  is 
stored  and  then  appropriately  weighted  by  the  field  radiated  from  St  of  the 
first  section  which  is  known  from  the  first  part.  The  resultant  information 
is  required  in  a  reciprocity  integral  formulation  which  then  formally  con¬ 
nects  the  appropriate  fields  associated  with  the  two  separate  parts  in  this 
hybrid  procedure  to  yield  the  desired  external  field  scattered  by  the  com¬ 
plex  interior  obstacle  in  the  original  configuration.  The  main  advantage  of 
this  hybrid  approach,  which  has  been  recently  shown  to  be  successful  by 
considering  several  representative  examples,  is  that  it  allows  one  to  sepa¬ 
rately  estimate  the  effects  on  the  overall  scattering  due  to  changes  in  the 
cavity  shape  for  a  given  type  of  complex  interior  obstacle,  and  due  to  differ¬ 
ent  complex  interior  obstacles  for  a  given  cavity  shape,  respectively.  This 
work  was  also  supported  in  part  by  General  Electric  of  Cincinnati,  Ohio. 
Additional  extensions  of  this  procedure  are  under  investigation. 

c.  Hybrid  Analysis  of  EM  Scattering  by  Antenna  Cavity  Con¬ 
figurations 

The  EM  scattering  by  antenna  cavity  shapes  is  being  studied  using  hybrid 
methods  under  JSEP  support.  Previous  work  under  JSEP  included  a  hy¬ 
brid  combination  of  asymptotic  high  frequency  and  modal  techniques  for 
the  analysis  of  the  scattering  from  two-dimensional  (2-D)  antenna  cavity 
shapes  for  which  the  modes  can  be  found  analytically  as  well  as  a  hybrid 
moment  method/asymptotic  method  analysis  for  2-D  cavities  with  domi¬ 
nant  mode  waveguide  fed  antennas.  Current  effort  is  focussed  on  hybrid 
methods  for  the  analysis  of  three  dimensional  (3-D)  antenna  cavity  config¬ 
urations. 
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This  past  .year,  the  analysis  of  the  scattering  from  a  3-D  rectangular 
antenna  cavity  backed  by  an  array  of  waveguide  fed  slots  was  completed. 
This  analysis  can  include  the  effect  of  dielectric  material  loading  on  both  the 
waveguide  as  well  as  the  cavity.  Additionally,  the  dielectric  loading  material 
may  be  different  for  the  two  different  regions.  This  cavity  configuration 
was  analyzed  in  two  different  ways,  using  a  full  moment  method  as  well  as 
using  an  analytical  approach  based  on  the  single  mode  approximation  for 
dominant  mode  waveguide  fed  antennas  that  allow  such  an  approximation. 
Both  methods  were  found  to  be  in  good  agreement. 

Work  is  currently  underway  for  the  analysis  of  generally  shaped  3-D 
antenna  cavities.  Using  the  method  of  moments,  the  mutual  coupling  and 
self  admittance  between  waveguide  fed  antennas  without  the  presence  of  the 
cavity  walls  has  been  determined.  Work  is  in  progress  to  find  the  coupling  of 
the  antennas  to  arbitrarily  shaped  cavity  walls  by  applying  high  frequency 
asymptotic  methods. 

d.  Hybrid  Analysis  of  EM  Scattering  from  Complex  Structures 

A  hybrid  combination  of  asymptotic  high  frequency  and  moment  method 
techniques  is  being  developed  to  analyze  the  EM  scattering  from  electri¬ 
cally  large  complex  structures.  Asymptotic  methods  are  used  on  objects 
which  are  electrically  large  and  able  to  be  locally  approximated  by  canoni¬ 
cal  shapes.  However,  at  lower  frequencies,  such  methods  are  not  accurate. 
On  the  other  hand,  the  moment  method  is  better  suited  for  analysis  of  ob¬ 
jects  which  are  electrically  small.  It  is  a  highly  accurate  method  but  very 
inefficient  especially  at  high  frequencies.  The  hybrid  method  is  designed 
to  be  a  bridge  between  these  two  methods.  Essentially  it  uses  asymptotics 
whenever  it  can  to  represent  just  a  few  moment  method  basis  functions  for 
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the  unknown  current  on  the  electrically  large  parts  of  the  structure,  and 
j.ict  a  few  conventional  moment  method  suhsectinnal  basis  functions  on  the 
electrically  small,  irregular  parts  of  the  geometry.  Thus,  the  advantage  of 
the  hybrid  method  lies  in  its  ability  to  analyze  an  object  which  is  electri¬ 
cally  large  but  arbitrary  without  having  to  use  an  extremely  large  number 
of  conventional  subsectional  or  even  entire  domain  basis  functions  over  the 
whole  radiating  object. 

The  configuration  being  presently  studied  is  a  finite  perfectly-conducting 
plate  of  relatively  arbitrary  shape  on  an  electrically  large  perfectly-conducting 
convex  cylinder.  Such  a  configuration  can  simulate,  for  example,  a  fin  on 
a  fuselage  of  an  aircraft  or  a  missile.  This  geometry  is  too  large  to  be  ana¬ 
lyzed  efficiently  by  the  moment  method;  yet,  no  asymptotic  solution  exists 
for  it.  For  this  geometry,  the  conventional  subsectional  basis  functions  will 
be  used  only  to  represent  the  unknown  current  on  the  fin  which  is  assumed 
to  be  electrically  small  while  the  cylinder  contribution  will  be  taken  into 
account  by  a  few  basis  functions  based  on  high  frequency  asymptotics  (or 
UTD). 

To  generalize  this  method  even  further,  we  will  consider  how  to  treat 
a  fin  that  is  not  electrically  small.  At  present,  we  contemplate  using  sub¬ 
sectional  basis  set  in  the  regions  where  the  current  is  rapidly  varying,  i.e. 
corners  and  edges,  and  asymptotics  elsewhere.  This  allows  many  more 
complex  objects  to  be  analyzed  by  such  a  hybrid  procedure. 
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VI.  ADAPTIVE  ARRAY  STUDIES 


Researchers: 

R.  T.  Compton,  Jr.,  Professor  (Phone:  614-292-5048) 

James  Ward,  Graduate  Research  Associate  (now  at  MIT  Lincoln  Lab) 

1.  Introduction 

During  the  past  year  we  have  done  research  under  the  JSEP  contract  in 
two  areas  related  to  arrays: 

1.  High  throughput  packet  radio  with  adaptive  arrays 

2.  High  resolution  angle  estimation 
We  discuss  these  areas  below. 

2.  High  Throughput  Packet  Radio  with  Adaptive  Ar¬ 
rays 

We  have  continued  our  research  on  adaptive  arrays  in  packet  radio  sys¬ 
tems  from  last  year.  We  have  obtained  exciting  results  showing  that  the 
throughput-delay  performance  of  random  access  packet  radio  systems  can 
be  substantially  improved  by  using  adaptive  arrays. 

At  the  end  of  last  year’s  JSEP  contract,  we  had  completed  an  initial 
study  of  the  performance  of  a  single-hop  slotted  ALOHA  system  using  a 
single-beam  adaptive  array  and  were  preparing  a  paper  on  that,  work.  (That 
paper  was  submitted  this  year  [1].)  We  had  also  begun  work  on  the  use  of 
a  multiple-beam  adaptive  array  (MBAA)  in  a  slotted  packet  radio  system. 
This  work  has  been  continued  during  the  past  year,  and  a  paper  describing 
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an  MBAA  in  a  slotted  system  and  deriving  its  performance  has  also  been 
submitted  [2]. 

In  most  packet  radio  protocols  studied  to  date,  it  is  assumed  that  a 
given  terminal  in  a  packet  radio  network  can  receive  only  one  packet  at  a 
time.  When  two  packets  collide,  it  is  assumed  that  neither  packet  is  received 
correctly.  In  some  cases,  a  capture  effect  [3, 4, 5, 6, 7]  may  exist  due  to,  for 
example,  differences  in  power  levels  of  the  incoming  packets.  A  capture 
effect  allows  one  packet  to  be  received  correctly  even  if  a  collision  occurs. 
But  even  with  capture,  it  is  still  assumed  that  a  terminal  can  receive  only 
one  packet  at  a  time. 

However,  this  limitation  can  easily  be  overcome  by  using  a  multiple* 
beam  adaptive  array  in  the  packet  system.  An  MBAA  is  an  adaptive 
antenna  system  in  which  several  array  output  signals  are  formed  simul¬ 
taneously  from  the  same  antenna  elements.  Each  output  signal  is  obtained 
by  using  a  different  set  of  adaptive  weights  to  combine  the  element  signals. 
In  this  way,  each  output  signal  has  a  different  pattern,  with  its  beam  and 
nulls  in  different  directions. 

By  using  an  MBAA  in  a  packet  system,  it  is  possible  to  receive  multiple 
packets  simultaneously  in  each  slot.  The  antenna  patterns  can  be  adap¬ 
tively  controlled  so  each  beam  of  the  MBAA  has  a  maximum  pointed  at  a 
different  incoming  packet  and  has  nulls  on  colliding  packets.  This  technique 
yields  enormous  increases  in  throughput  and  reductions  in  delay. 

The  performance  of  a  single-beam  adaptive  array  in  a  packet  system  is 
similar  to  that  of  Carrier  Sense  Multiple  Access1  (C'SMA)  [8] .  But  when 
a  multiple-beam  adaptive  array  is  used,  a  much  more  substantial  improve¬ 
ment  is  obtained.  One  can  achieve  average  throughputs  of  2  to  3  packets 

'However,  an  adaptive  array  has  the  advantage  over  CSMA  that  with  an  adaptive  array 
it  is  not  necessary  that  all  terminals  in  the  network  be  able  to  hear  each  other. 


44 


per  slot,  for  example,  with  only  modest  MBAA  capabilities. 

The  Mse  of  r.n  MB  A  A  to  achieve  high  throughput  and  low  delay  in  a 
packet  system  may  be  viewed  as  an  alternative  to  spread-spectrum  tech¬ 
niques,  which  have  also  been  considered  as  a  means  of  receiving  more  than 
one  packet  in  the  same  slot  [9,10].  The  use  of  spread  spectrum  also  allows 
multiple  packets  to  be  received,  either  by  separating  packets  on  the  basis 
of  arrival  times  (with  delay  capture,  which  separates  packets  with  the  same 
code  [11])  or  by  using  orthogonal  codes  on  different  packets  (with  Code 
Division  Multiple  Access  (CDMA)  [9,10]).  However,  spread  spectrum  tech¬ 
niques  have  the  disadvantage  that  they  require  much  wider  bandwidths 
than  those  needed  to  support  the  data  modulation.  An  MBAA  system 
requires  only  the  data  bandwidth. 

Figure  10  shows  a  typical  performance  curve  for  a  packet  radio  receiver 
with  an  MBAA.  It  shows  the  average  conditional  throughput  S(j)  as  a 
function  of  the  number  of  blocked  users  j  for  a  network  with  50  terminals 
and  an  MBAA  with  K  =  1  to  6  beams.  These  curves  were  computed 
for  a  new  packet  transmission  probability  Pn  =  0.02,  a  blocked  packet 
retransmission  probability  PT  =  0.2,  and  for  an  array  with  10  elements 
that  can  resolve  packets  5°  apart.  It  may  be  seen  how  the  conditional 
throughput  increases  significantly  when  the  adaptive  array  is  included  and 
as  the  number  of  beams  is  increased.  Conditional  throughputs  of  nearly 
4  packets  per  slot  are  possible  with  a  6-beam  adaptive  array.  Adding  an 
MBAA  to  a  packet  system  also  improves  the  stability  of  the  system. 

During  this  year,  we  have  also  begun  studying  the  use  <>f  adaptive  ar¬ 
rays  in  unslotted  ALOHA  systems.  An  unslotted  system  is  more  difficult 
to  analyze  than  a  slotted  system.  Also,  an  MBAA  in  an  unslotted  system 
requires  quite  different  techniques  for  timing  acquisition  and  array  adap- 
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S(j),  the  throughput  at  state  j 


Figure  10:  Average  conditional  throughput  S(j )  and  traffic  load  Sm{j) 
for  a  network  with  50  terminals  and  an  MBA  A  with  K  beams.  (New 
packet  transmission  probability  Pn  =  0.02,  backlogged  packet  transmission 
probability  Pr  =  0.2,  8  nulls,  63  bit  PN  acquisition  code.) 
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tatioii  than  one  in  a  slotted  system,  because  of  the  asynchronous  nature 
of  the  unslotted  case.  Interfering  packets  may  arrive  at  any  time  during  a 
received  packet,  so  the  array  must  be  able  to  adapt  throughout  an  acquired 
packet.  (In  the  slotted  case,  adaptation  is  completed  during  the  preamble 
and  the  interfering  packets  are  nulled  prior  to  the  message  portion  of  the 
packet.) 

We  have  developed  a  promising  method  for  using  an  MBAA  in  an  un¬ 
slotted  system.  For  this  method  we  have  determined  the  throughput  of  an 
unslotted  ALOHA  receiver  using  an  MBAA  as  a  function  of  the  number 
of  beams,  the  number  of  degrees  of  freedom  and  the  angular  resolution  of 
the  array.  Our  analysis  is  based  on  a  method  similar  to  that  of  Tobagi 
[12]  as  applied  to  a  spread  spectrum  multiple  access  network  [10].  We  have 
also  done  simulations  to  verify  the  analysis.  Our  results  show  that  the 
throughput  of  an  unslotted  system  can  also  be  significantly  improved  with 
an  MBAA.  Moreover,  an  MBAA  brings  the  performance  of  an  unslotted 
system  close  to  that  of  a  slotted  system. 

Fig.  11  shows  the  throughput  for  a  network  of  M  —  20  terminals  when 
the  packet  receiver  uses  an  adaptive  array  with  a  single  beam.  Curves 
are  shown  for  different  numbers  of  nulls  N  in  the  array.  The  N  =  0  case 
corresponds  to  standard  unslotted  ALOHA  (with  zero  capture).  These 
curves  are  for  the  ideal  case  when  the  array  can  resolve  signals  arbitrarily 
close  in  angle.  The  curve  for  N  ~  19  corresponds  to  perfect  capture. 
For  all  traffic  loads  G,  the  throughput  increases  as  N  is  increased.  The 
improvement  in  throughput  with  additional  nulls  may  lie  seen. 

Fig.  12  shows  typical  curves  of  the  throughput  of  an  unslotted  system 
with  A/  =  20  terminals  using  an  MBAA  with  N  8  nulls.  The  parameter 
K  in  the  figure  is  the  number  of  beams  formed  by  the  MBAA.  Curves  are 
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shown  for  K  =  1.3,5,  7  and  9.  In  Fig.  12a,  the  array  can  resolve  packets 
with  arbitrarily  small  angular  separation  ( 0r  =  0°),  but  in  Fig.  12b  the 
array  can  resolve  packets  only  above  a  ffT  =  5°  separation.  Fig.  12b  shows 
theoretical  results  in  the  solid  curves  and  simulation  results  in  the  dashed 
curves.  Note  that  this  system  achieves  an  average  throughput  above  2.5 
packets  per  packet  length,  and  it  does  this  in  an  unslotted  system! 

Finally,  Fig.  13  shows  a  comparison  between  a  slotted  and  an  unslotted 
system  using  an  array  with  8  nulls  in  ::  network  of  20  terminals. 

3.  High  Resolution  Angle  Estimation  Techniques  in 
Packet  Radio 

Our  work  has  shown  that  the  performance  of  packet  radio  networks  can  be 
significantly  improved  by  taking  advantage  of  the  directional  properties  of 
arrays.  This  observation  has  motivated  us  to  consider  other  ways  of  using 
active  antenna  processing  in  packet  systems. 

It  appears  that  an  even  greater  improvement  in  packet  radio  throughput 
and  delay  could  be  achieved  by  using  angle  estimation  techniques  in  such 
systems.  For  example,  consider  a  slotted  packet  radio  system  in  which 
an  array  with  several  elements  L  used  at  a  packet  receiver.  Suppose  that 
instead  of  receiving  packets  in  real  time  during  each  slot,  we  instead  sample 
the  output  of  each  array  element  during  the  slot  interval  and  store  the 
samples  in  computer  memory.  Then,  at  the  end  of  each  slot,  after  the  signals 
have  been  received,  we  combine  the  samples  from  the  array  elements  with 
several  different  sets  of  array  weights.  Each  set  of  weights  will  be  chosen  In 
receive  one  packet  and  null  the  others.  In  order  to  do  this,  il  is  necessary 
to  do  two  things  at  the  end  of  each  slot: 

1.  Estimate  the  number  of  packets  J  that  arrived  during  the  slot,  and 
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Throughput  S 


Figure  11:  Throughput  for  a  single-beam  adaptive  array.  M  =  20,  0T  =  0° 
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(a) 


0>) 


Figure  12:  Throughput  for  an  MBAA  with  N  ~  8  nulls  and  M  —  20 
terminals,  (a)  ffr  —  0°.  (h)  9r  -  ,r)°.  Solid  line:  theory;  Dashed  line: 

simulation) 
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Figure  13:  Comparison  of  unslotted  and  slotted  throughput  for  a  10-ele 
ment  MBAA  as  the  number  of  beams  is  varied.  M  =  20,  N  =  8,  8r  —  0° 
(Solid:  slotted;  dashed:  unslotted-exponontial  packet  lengths;  dashdot:  un 
slotted-fixed  packet  lengths) 
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2.  Estimate  the  arrival  angles  Qj,  1  <3<J,  for  each  of  the  J  packets 
in  the  slot. 

When  J  and  the  Oj  have  been  estimated,  we  then  calculate  J  sets  of 
array  weights.  Each  set  of  array  weights  will  be  chosen  to  form  a  pattern 
with  its  maximum  toward  one  of  the  incident  packets  and  nulls  on  all  the 
others.  These  array  weights  can  then  be  applied  to  the  signal  samples 
stored  in  RAM  to  create  J  different  array  output  signals,  each  containing 
only  one  incident  packet.  In  this  way,  all  the  incident  packets  in  each 
slot  can  be  received  without  distortion.  The  cost  required  to  achieve  this 
advantage  is  the  fixed  delay  involved  in  storing  the  samples  and  weighting 
them  after  the  end  of  the  slot.  However,  this  fixed  delay  is  small  compared 
to  the  retransmission  delay  encountered  in  a  busy  network  that  is  trying  to 
receive  packets  in  real  time.  The  ultimate  performance  of  such  a  “delay- 
and-weight”  system  depends  on  the  number  of  array  elements,  the  packet 
signal-to-noise  ratios,  the  angular  separation  between  packets,  and  so  forth. 

During  the  past  year  of  the  JSEP  contract,  we  have  begun  to  study  this 
approach.  We  started  by  considering  the  ang  estimation  techniques  that 
would  be  used  to  implement  such  a  system.  To  date,  we  have  submitted 
two  papers  for  publication  based  on  this  work. 

The  first  paper  [13]  considers  Maximum  Likelihood  (ML)  estimation  of 
packet  arrival  angles  for  the  case  where  each  packet  contains  a  known  PN 
(pseudonoise)  code  at  the  beginning  of  the  packet.  Such  an  acquisition 
code  was  used  for  our  MBA  A  technique  and  we  were  interested  in  whether 
including  such  an  acquisition  code  improves  the  performance  of  an  ML 
estimator. 

In  a  second  paper  [14],  we  have  considered  the  use  of  the  ESPRIT 
algorithm  [15[  and  shown  how*  it  may  be  used  in  a  cross-polarized  arrays  to 
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estimate  both  the  arrival  angles  and  the  polarizations  of  signals. 

We  plan  to  continue  work  in  this  area  under  the  follow-on  JSEP  con¬ 
tract. 
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